Anisotropic spin-dynamics in the Kondo semiconductor CeRu 2 Ali 
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Spin dynamics in the new Kondo insulator compound CeRu2Alio has been studied using unpo- 
larized and polarized neutron scattering on single crystals. In the unconventional ordered phase 
forming below To = 27.3 K, two excitation branches are observed with significant intensities, the 
lower one of which has a gap of 4.8 ± 0.3 meV and a pronounced dispersion up to w 8.5 meV. Com- 
parison with RPA magnon calculations assuming crystal-field and anisotropic exchange couplings 
captures major aspects of the data, but leaves unexplained discrepancies, pointing to a key role of 
direction-specific hybridization between 4/ and conduction band states in this compound. 

PACS numbers: 71.27.+a, 75.20.Hr, 75.30.Ds, 75.30.Gw, 75.30.Mb, 78.70.Nx 
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The CeM 2 Ali (M = Fe, Ru, Os) compounds form 
a new family of Ce-based intermetallic materials with 
fascinating, but hitherto elusive, magnetic and trans- 
port properties. Below room temperature, they show 
evidence of a Kondo-insulator regime, with an increase 
in the electrical resistivity on cooling ascribed to the 
opening of a narrow "hybridization gap" in the elec- 
tronic density of states [U, Hj] ■ In the standard approach 
Q, this mechanism should ultimately lead to a non- 
magnetic, many-body singlet ground state for T — >• 0, 
as was observed experimentally for the vast majority 
of Kondo-insulator compounds known to date. In con- 
trast, CeRt^Alio and CeOs2Alio order magnetically be- 
low T = 27.3 K and 28.7 K, respectively Q. Their struc- 
ture is antiferromagnetic (AF) with the simple wavevec- 
tor kAF = (0,1,0) However, there is strong ex- 

perimental evidence that this ordering cannot be ex- 
plained by conventional Ruderman-Kittel Kasuya-Yosida 
(RKKY) exchange alone: To seems unrealistically high in 
view of the large Ce-Ce interatomic distances (5.26 A), 
of the weak ordered antiferromagnetic moment ([aaf = 
0.32(4)-0.4/xb 0-01 for M = Ru) derived from neu- 
tron diffraction measurements, and of the much lower 
Neel temperatures found in other TRU2AI10 compounds 
(T N = 16.5 K in GdRu 2 Al 10 [8]). It was also reported 
that T increases with the application of pressure [l|, 
contrary to the general trend in Ce Kondo compounds. 
This unique situation has attracted considerable inter- 
est because it seems to challenge widely accepted views 
on Kondo insulators. Various interpretations have been 
proposed in terms of (i) a charge density wave associ- 
ated with an energy gap opening preferentially along the 
b direction 0, [l(| , (it) a spin-Peierls state due to the for- 
mation of spin-singlet pairs [lll4l3j . or (Hi) a resonating- 
valence-bond state [14j |. Quite remarkably, despite the 



large anisotropy of the paramagnetic susceptibility with 
Xa S> Xc 3> Xbt the ordered AF moments align along 
the c direction @, @|. In Refs. [H, EH, this discrep- 
ancy was suggested to arise from conduction-electron- 
/-electron (c—f) hybridization occurring predominantly 
along a, and suppressing \a accordingly through the for- 
mation of a (Kondo) spin singlet. A detailed study of 
the spin dynamics is of primary importance to sort out 
this problem. Inelastic neutron scattering (INS) experi- 
ments performed previously on powder samples have ev- 
idenced the opening of a large spin-gap in the ordered 
state, with a broad excitation centered at Ago = 8 meV 



and 11 meV in CeRu2Ali [4] and CeOs 2 Ali [17], respec- 
tively. However, mode dispersion and anisotropy were 
obscured by powder averaging, and it could not be de- 
cided whether the observed magnetic signal arose from 
dispersive magnon branches with an anisotropy gap or, 
e.g., from singlet-triplet transitions with sizable disper- 
sion and/or damping. The possibility of a lattice con- 
tribution could also not be ruled out. In this Letter, 
we report unpolarized and polarized INS experiments 
performed on single-crystal CeRu2Alio- The spectra re- 
veal well-defined dispersive excitations with a gap of 4.8 
meV at the AF zone center. They exhibit a remarkable 
anisotropy which does not correspond to a standard pre- 
cession of spin wave modes. Overall agreement with the 
experimental results can be achieved phenomenologically 
in a RPA model by assuming a strongly anisotropic bi- 
linear exchange interaction J c ^> J a ,J h . However, re- 
maining inconsistencies are thought to reflect anisotropic 
hybridization effects, whose role was suspected from pre- 
vious studies (lif . 

Thirteen single crystals of CeRu2Alio (orthorhombic, 
Cmcm space group, No. 63) with dimensions comprised 
between 1 and 4 mm, for a total mass of about 500 mg, 



2 



800 



600 



(a) 


Q = (1,0,1) 


< b > t=iok 




* 

A 


Q= (0,3,0) M 




L / \ 7=3.2K 


i / T ■ 












5| (3, 0, 0) ^ 








Subtracted , / \/ 






background \ J 





1200 



1000, 



10 



800 o 



02468 10 02468 10 12 
Energy (meV) 

FIG. 1. (Color online) Energy scans measured on 2T (a) and 
IN8 (b) at k f = 2.662 A" 1 for different AF zone centers: (a) 
Q = (1, 0, 1) at T = 3.2 K (closed circles) and 38 K (open cir- 
cles). A steep, temperature independent, background (dashed 
line) was estimated by assuming the magnetic signal at 38 
K to be Q independent (verified for other Q vectors), then 
subtracted from the measured data, (b) Q = (3, 0, 0) (open 
circles) and (0, 3, 0) (closed squares) at T — 10 K. In (a) 
and (b), solid lines represent fits to the data using Gaussian 
(inelastic) and Lorentzian (quasielastic) line shapes. 



were grown by an Al-flux method, and co-aligned with 
their b axes vertical on an Al sample holder. An effec- 
tive mosaicity of about 3 degrees was estimated from 
the neutron rocking curves, which was sufficient for the 
present experiment. Excitation spectra were measured in 
the (o*,c*) scattering plane, first using unpolarized neu- 
trons on the 2T triple-axis spectrometer at LLB-Orphec 
(Saclay), then with linear polarization analysis on IN20 
at the ILL (Grenoble). Finally, the crystals were reori- 
ented with the c axis vertical on a lighter sample holder, 
and measured with unpolarized neutrons on IN8 (ILL) 
in a 6 T cryomagnet. Spectra were recorded at fixed fi- 
nal energy, Ef = 14.7 meV, using a pyrolytic graphite, 
PG002 (2T) or Silll (IN8) monochromator and a PG002 
analyzer, with a PG filter placed on the scattered beam, 
or (polarized neutrons on IN20) a Heussler monochroma- 
tor and analyzer. 

Constant-Q scans have been performed using unpolar- 
ized neutrons for momentum transfers lying in the (a*, 
c*) (2T) and (a*, b*) (IN8) planes. Representative spec- 
tra are presented in Fig. Q] For T = 3.2 K (2T) or 10 
K (IN8), one or two distinct modes are visible depending 
on the Q vector. The dispersion is significant, with a gap 
of 4.8 ± 0.3 meV at the AF zone centers (Fig. [2}. Near 
the zone boundary, the excitations reach 8.5 ± 0.3 meV, 
with a flat region corresponding to the peak observed just 
above 8 meV in the previous powder experiments Q • In- 
tensity maps for three particular directions, Q = (1, 0, 1), 
(h, 0, 1) and (h,3, 0), are presented in Fig. [3] The exis- 
tence of (at least) two modes is best evidenced in scans 
at q = k A F, e.g. for Q = t 02 o + k A F = (0, 3, 0). On the 
other hand, the lower branch shows no detectable inten- 
sity at the AF Q vector (1, 0, 2) (of the form (h, 0, /) with 
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FIG. 2. (Color online) Dispersion of the magnetic excitations 
in CeRu2Alio. Closed circles: results from the unpolarized 
neutron experiments; red and blue symbols: guides to the 
eye; grey lines: RPA calculations (see text). 
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FIG. 3. (Color online) Left: intensity maps derived from the 
spectra measured on (a, b): 2T at T = 3.2 K and (c): IN8 at 
T = 10 K for three directions in reciprocal space. Different 
color scales are used to reflect the different counting rates on 
the two spectrometers. Right: RPA calculations (see text). 



I even), as shown in Fig. [3^a) . Another important obser- 
vation is that the magnetic intensity of the lower branch 
is strongly suppressed for scattering vectors whose orien- 
tation is close to the a* axis, such as Q = (3, 0, 0) as 
compared to (0, 3, 0) (Fig.Q]). This suggests that dynam- 
ical correlations (m^mfy and (m°m^) between moment 



3 



300 



200 



100 



(1,0, 3) 
(3, 0, 1) 




8 10 2 
Energy (meV) 

10 




FIG. 4. (Color online) Polarization analysis on IN20 (kf — 
2.662 A -1 ). Upper frames : energy scans measured at T — 2 
and 45 K for different AF zone centers. Left: Q = (1,0, 1); 
grey symbols denote the NSF intensity measured with Po || x, 
and other symbols SF intensities for P || x, y or z. Right: 
Q = (1, 0, 3) and (3, 0, 1) (red squares and blue diamonds, re- 
spectively); the plot shows SF intensities for P || x. Full 
(dashed) lines represent intensities calculated using Gaus- 
sian (Lorentzian) spectral functions. Lower frames: intensity 
maps along the (h, 0, 1) direction for the two transverse com- 
ponents {mfrrij} (left) and (mfm^} (right) of the magnetic 
correlations using the exchange parameters listed in Table U 



components perpendicular to the a axis are weak, and 
(mfrrij} correlations dominate the magnetic response. 
When temperature increases to T — 38 K > To, the in- 
elastic magnetic peak at 4.8 ± 0.3 meV is suppressed and 
replaced by a sloping intensity at low energy. The latter 
signal shows no pronounced Q dependence [apart from 
the appearance of a strong extra background near Q = 
(1, 0, 1), see Fig.QJa)], and is thus ascribed to quasielas- 
tic (QE) fluctuations. Spectra along (h, 3, 0), (0, 2 + k, 0), 
and (h, 3 — ft., 0) were also measured at the base tempera- 
ture in an applied field (H || c) of 5 T, above the moment 
reorientation transition from mAF || c to mAF || b, known 
to occur at 77* » 4 T 3, 3j. No sizable change was 



observed with respect to the H = data. 

Neutron polarization analysis provides further insight 
into the anisotropy of the magnetic response. Fig. [4] 
(upper frames) shows intensities measured in the spin- 
flip (SF) and non-spin-flip (NSF) channels at different 
scattering vectors. Let us first consider the results for 
Q = (1,0,1). One sees that the NSF signal (measured 
with the incident polarization P || Q) is featureless and 
temperature independent, confirming the magnetic origin 
of both the 4.8-meV peak below 10 K and the QE signal 



at T > Tb found in the unpolarized neutron experiments. 
The SF intensities measured for different directions of in- 
cident polarization P || x,y or z [13] are found to fulfill 



I x ~ I z and I y ~ 0. Using the standard expressions [21 1 



if oc M 00 (q) sin 2 (a) + M bb (q) + M cc (q) cos 2 (a)] (la) 



if oc M 66 (q) 

if cc M aa (q) sin 2 (a) + M cc (q) cos 2 (a), 



(lb) 
(1c) 



where M pp is the dynamic structure factor associated 
with pair correlations of the moment component m p 
(p = {a, b, c}), and a the angle between Q and the a* 
axis, one comes to the conclusion that correlations of 
the nib components must vanish to the precision of the 
present measurement. For Q = (1,0, 1), a is very close 
to 45° since the lattice parameters a and c are nearly 
equal. In contrast, the scattering vectors Q = (1,0,3) 
and (3, 0, 1) correspond to the same reduced q vector 
(AF zone center) and nearly equal values of the dipole 
magnetic form factor, but their a angles are quite dif- 
ferent (71.6° and 18.4°, respectively). From Fig. [4] the 
ratio of the magnetic excitation intensities 7^ for those 
two spectra is about 3.25, which implies that correlations 
of a components dominate. Assuming Mbb to be strictly 
zero, and solving Eqns . fTal and [Tel one gets M aa /M cc 5. 
In a magnon picture, such a difference can be understood 
by noting that M aa and M cc correspond, respectively, to 
transverse and longitudinal excitation modes of the AF 
magnetic structure. On the other hand, the strong differ- 
ence between the transverse components along a and b is 
quite remarkable and requires a very unusual anisotropy 
to exist in this material. 

To analyze this magnetic response, we have per- 
formed calculations assuming bilinear exchange interac- 
tions, Hi j = ^aJ^SfS", between near-neighbor 
Ce sites. Both a standard spin wave model, and random- 
phase approximation (RPA) calculations were investi- 
gated. In the following, we will focus on the second 
approach, which can treat anisotropy effects in a more 
realistic way. The crystal-field (CF) parameters for the 
Ce 3+ J = 5/2 ground state, (B%, 77f, B%, 77f, Bf) = 
(-1.326, -29.236, +1.013,-1.747, -5.317) K, choosing c 
as the quantization axis, were taken from Strigari's work 
[22l |. and correspond to a sequence of three doublets at 0, 
354 K, and 535 K. The resulting single-ion anisotropy has 
an easy a axis, as required by the magnetic susceptibility 
measured in the paramagnetic regime. Therefore, in this 
simple picture, one has to assume that J c is much larger 
than J a and J b to ensure that the AF ordered moments 
properly align along the c axis. 

Fair overall agreement can be obtained between the 
calculations and the experimental excitation spectra be- 
low To, using the set of exchange constants listed in Ta- 
ble UJ As can be seen in Fig. [31 the observation of two 
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TABLE I. Anisotropic exchange parameters (in units of K) 
used in the RPA calculation. Atomic positions (xi,yi, Zi 
1 : (0,2/, i);2 : (1/2, U-y, |);3 : (1/2, | -y, |),4 : (0, - 
with i/ = 1.1239(3) HI] 
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branches with significant spectral weight (from a total of 
4) , as well as the general Q dependence of their intensities 
along different symmetry directions, or the anisotropy of 
the correlations (Fig. HI lower frames) can be accounted 
for. Furthermore, salient features of the experimental 
data are well reproduced in the calculations, such as the 
vanishing of the lower dispersive branch near the (1,0,2) 
AF zone center (upper frames in Fig. [3j) , or the significant 
intensity exhibited by the upper branch near Q = (0, 3, 0) 
[lower frames, in accordance with Fig.[T](b)], in contrast 
to, e.g., Q = (1,0,1) [upper frames and Fig. Ufa)]. On 
the other hand, notable quantitative differences exist: the 
initial slopes of the dispersions are much steeper than 
predicted by the calculation, and the calculated energy 
of the upper mode is too high. We believe that this dis- 
crepancy results from the unrealistically large J c value 
required to keep the ordered moments aligned along the 
c axis despite the strong single-ion anisotropy favoring 
the a axis. Simulations done in the simpler Holstein- 
Primakoff spin-wave approximation indeed showed that 
the agreement improves if one reduces this single-ion 
anisotropy and, correspondingly, the anisotropic compo- 
nent of the exchange tensor. Recent simulations per- 
formed in a mean-field, two-sublattice, model [24[ fur- 
ther indicate that anisotropic exchange parameters large 
enough to overcome the single-ion a-axis anisotropy in- 
evitably result in a large ordered moment, contrary to the 
experimental observation of \iaf — 0.32(4)-0.42(1)^b. 
This could raises the question of whether the CF model 
of Ref . 22| used in the present calculations overestimate 
the single-ion anisotropy. Meanwhile, there is growin g ex - 
perimental evidence, as discussed in recent papers [jl Il6l| , 
that direction-selective hybridization of 4/ orbitals with 
conduction band states plays a key role in the pecu- 
liar magnetism of the CeT2Alio compounds. This has 
been proposed to explain the anomalous magnitude of 
the single- ion anisotropy in CeRu2Alio, as compared to 
that of NdT2Alio, as well as the lack of a sizable anomaly 
(Al/l < 10~ 6 ) in the longitudinal magnetostriction at the 



critical field H^ c «4T where the AF moment direction 
reorients from c to b (possibly related to the intrigu- 
ing lack of field dependence of the magnetic excitation 
spectra found in the present measurements) . It has been 
argued [23] that, owing to specifics of the YbFe2Alio-type 
crystal structure, f—p hybridization takes place predom- 



inantly within the (a, c) plane, especially with the Al(2) 
atoms located in the a direction with respect to the Ce 
site. This hybridization could result in a suppression of 
the magnetic components along a, thereby favoring the 
alignment of the ordered AF moments along c. Such a 
picture provides an appealing physical basis for the re- 
duction of the single-ion anisotropy hypothesized in the 
above discussion. In the case of CeOs2Ali , it has been 
argued [HI that the gap in the magnetic excitation spec- 
trum, associated with the formation of a singlet state, 
starts to develop below the temperature of the maximum 
in the magnetic susceptibility \a (T) , well above the onset 
of the AF order. Such effects are clearly beyond the scope 
of the simple magnon model presented above, which basi- 
cally treats the spin gap as an anisotropy gap, but should 
be included in a more realistic treatment. 

In conclusion, the present study provides detailed in- 
sight into the spin dynamics of CeRu2Alio, and empha- 
sizes the most peculiar anisotropy of the magnetic corre- 
lations occurring in the AF ordered state. The results 
could be partly accounted for using a magnon model 
treated in the RPA approximation. However quantita- 
tive discrepancies suggest that this picture should be re- 
garded as phenomenological, and support the idea that 
anisotropic c-f hybridization plays a key role in this ma- 
terial. Proper theoretical consideration of such effects 
should open the way to a unifying view of static and dy- 
namic aspects of magnetism in this family of compounds. 
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